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Abstract 

In the present paper, we develop geometric analysis techniques on Cay- 
ley graphs of finitely generated abelian groups to study the polynomial 
growth harmonic functions. We provide a geometric analysis proof of the 
classical Heilbronn theorem [22] and the recent Nayar theorem [41] on 
polynomial growth harmonic functions on lattices Z" that does not use a 
representation formula for harmonic functions. In the abelian group case, 
by Yau's gradient estimate we actually give a simplified proof of the more 
general polynomial growth harmonic function theorem of Alexopoulos [3] 
for groups of polynomial volume growth. We also calculate the precise di- 
mension of the space of polynomial growth harmonic functions on finitely 
generated abelian groups. While the Cayley graph not only depends on 
the abelian group, but also on the choice of a generating set, we find that 
this dimension depends only on the group itself. 



1 Introduction 

Classically, in 1948 Heilbronn [35] proved the polynomial growth harmonic 
function theorem on the lattice Z" that polynomial growth discrete harmonic 
functions are polynomials, and calculated the dimension of the space of poly- 
nomial growth discrete harmonic functions. Recently, in 2009 Nayar [U] gave 
another proof of this theorem by the probabilistic method. Their proofs all 
depend on the representation formula for discrete harmonic functions. In this 
paper, we give a geometric analysis proof that can yield more general results. 

The study of harmonic polynomials on W 1 is classical, and the precise di- 
mension calculation of harmonic polynomials can be found in |33j . In 1975, Yau 
|5Qj proved the Liouville theorem for harmonic functions on Riemannian man- 
ifolds with nonnegative Ricci curvature. Then Cheng- Yau [5] used Bochner's 
technique to derive the gradient estimate for positive harmonic functions called 
Yau's gradient estimate which implies that sublinear growth harmonic functions 
on these manifolds are constant. Then Yau [5T][S2] conjectured that the space of 
polynomial growth harmonic functions with growth rate less than or equal to d 
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on Riemannian manifolds with nonnegative Ricci is of finite dimension. Li- Tarn 
[33] and Donnelly-Fefferman [TBI independently solved the conjecture for mani- 
folds of dimension two. Then Colding-Minicozzi jTQl [TTJ |T2] gave the affirmative 
answer by using the volume doubling property and the Poincare inequality for 
arbitrary dimension. The simplified argument by the mean value inequality can 
be found in |32[ 113) where the dimension estimate is asymptotically optimal. 
This inspired many generalizations on manifolds [3H1 [3FJ ESI 133 |3S1 13 I2Z1 135 
and on singular spaces [TOJ [2E1 [Ml US]- In this paper, we give the pre- 
cise dimension calculation of polynomial growth harmonic functions on finitely 
generated abelian groups. 

In another direction, Avellaneda-Lin [5] first proved the polynomial growth 
harmonic function theorem for elliptic differential operators with periodic coef- 
ficients in R™ (see [301 HI 121 SI 13H1 131] for more generalizations). Alexopoulos 
[3] proved a more general polynomial growth harmonic function theorem for 
groups of polynomial volume growth which contained our case. By a famous 
theorem of Gromov [20] , every finitely generated group G of polynomial growth 
is virtually nilpotent (i.e. it has a nilpotent subgroup H of finite index). For 
some torsion-free subgroup H ' of H, it can be embedded as a lattice in a simply 
connected nilpotent Lie group N. By considering the exponential coordinate of 
N, N is identified with R™ (see [42 ]). Alexopoulos proved that every polyno- 
mial growth harmonic function on G is the restriction to H' of some polynomial 
on R™. He used the homogenization theory and Krylov-Safonov's argument to 
compare the discrete heat kernel on H' and the continuous heat kernel for some 
sub-Laplacian on N. Then he adapted some ideas of Avellaneda-Lin [5] to prove 
a Taylor formula for harmonic functions on G. Instead of doing that, we shall 
prove Yau's gradient estimate on abelian groups to simplified the argument in 
this special case. 

We develop the geometric analysis techniques on Cayley graphs of finitely 
generated abelian groups. Firstly, we prove Yau's gradient estimate (see The- 
orem |TTT|) for positive discrete harmonic functions. Note that Kleiner [55] ob- 
tained the Poincare inequality on Cayley graphs of finitely generated (not nec- 
essarily abelian) groups (see also [45 J. For the abelian case, combining it with 
the natural volume doubling property, we obtain the uniform Poincare inequal- 
ity (see Lemma l3~Tj) from which the mean value inequality follows by the Moser 
iteration. In addition, for an abelian group, Bochner's formula for discrete 
harmonic functions can be easily calculated, i.e. |Vu| 2 is subharmonic for the 
discrete harmonic function u which is essentially due to Chung- Yau [9] and Lin- 
Yau [40] for Ricci fiat graphs. Then these results together imply Yau's gradient 
estimate on Cayley graphs of finitely generated abelian groups as in the case of 
Riemannian manifolds with nonnegative Ricci curvature. By the fundamental 
theorem of finitely generated abelian group (see [271 [33]), any finitely generated 
abelian group G is isomorphic to the direct sum Z m ©i=i^g< , where m, I G IN, 
qi = p"* for some prime number pi and ai 6 IN. 

Theorem 1.1 (Yau's gradient estimate). Let (G,S) be the Cayley graph of a 
finitely generated abelian group with symmetric generating set S (i.e. S = —S), 
and G = Z m @®\—i^q i - Then there exist constants C\ and Ci depending only 
on m and S, such that for any x G G, R > 1 and any positive discrete harmonic 
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function f on Bc 1 r{x) 1 we have 

|V/|(aO < ^f(x). (1.1) 

Secondly, from Yau's gradient estimate, we know that |V/| < Gi? d_1 , on 
Br for R > Rq if the discrete harmonic function / satisfies |/| < CR d on Br 
for R> Ri. That is, the growth order decreases when we take derivatives. 

This is the key to an induction argument to give a geometric analysis proof 
of Heilbronn's theorem. This scheme will, in fact, work for all abelian groups. 
Let us denote the set of polynomial growth harmonic functions of growth rate 
less than or equal to d on the Cayley graph (G,S) by H d (G,S) := {u : G — ► 
R | L s u = 0, |u|(x) < C{d s {p,x) + l) d }, where L s is the Laplacian operator 
on (G,S), d s is the distance function to some fixed p G G. For the finitely 
generated abelian group G — G\ © Gi = Z m ©| =1 Z gi , we define the natural 
projection 

7r Gl : G — > Gi, 

X >-» 7T Gl (x) = Xl, 

for x = x\ +X2, Xi G Gi and X2 G G2. It is easy to see that ttg 1 S is a generating 
set of Gi if S is a generating set of G. 

Theorem 1.2 (Generalized Heilbronn's theorem). Let (G,S) be the Cayley 
graph of a finitely generated abelian group, G — G\ © G2 — Z m 0©' =1 Z 9i . 
Then 

H d {G,S)=H d {G 1 ^ Gl S), 

moreover any f G H d (G,S) is a polynomial when it is restricted to G\ = W l , 
and it is constant on G2 i-e. f{x + w) = /(x), for \/x £ G, Vw G G2. 

Nayar |41j proved a strong version of Heilbronn's theorem. We denote by 
HM d (G,S) := Span{w : G -> R | L s u = 0,u(x) > -C(d s (p,x) + l) d } the 
linear span of one-sided bounded polynomial growth harmonic functions. It is 
trivial that H d (G,S) C HM d (G, S). By the Harnack inequality (03, we give 
a geometric analysis proof of Nayar's theorem. 

Theorem 1.3 (Generalized Nayar's theorem). Let (G,S) be the Cayley graph 
of a finitely generated abelian group, G — G\ © G2 = Z m ffi' = i^ 9i • Then 

HM d (G, S) = H d (G, S) = H d (Gi , ttq 1 S), 

moreover any f G H d (G,S) is a polynomial when it is restricted to G\ = Z m 
and it is constant on Gi- 

Thirdly, for discrete harmonic polynomials we calculate the precise dimen- 
sion by linear algebra. In fact, instead of using the technical lemma (see [T91I23] ) 
from difference equations as Heilbronn [35] did in the lattice case, we apply the 
dimension comparison argument with harmonic polynomials in R™. Conversely, 
our argument provides a proof of this difference equation lemma. Moreover, 
we can calculate the dimension of the space of polynomial growth harmonic 
functions on Cayley graphs of finitely generated abelian groups. It is surprising 
that the dimension of polynomial growth harmonic functions does not depend 
on the generating set S for the abelian group G. Actually, the graph structures 
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of two Cayley graphs of the abelian group G with two generating set S\ , S2 can 
be quite different. While the Laplacian operator depends on the generating set, 
the dimension of polynomial growth harmonic functions does not. We denote 
by HP k (W n ) the space of harmonic polynomials on E, m with degree less than 
or equal to k, k G M U {0}. It is well known (see [33 ) that 

m + k - l\ (m + k - 2\ 
k ) + { k-1 J' 

Theorem 1.4 (Dimension calculation). Let (G,S) be the Cayley graph of a 
finitely generated abelian group, G = G\ © G2 — ^ m @ ®' =1 Z gj . Then 

dim HM d (G, S) = dim H d (G, S) = dimHP [d] (R m ). 

Since we do not use the representation formula, this method can be ap- 
plied in more general settings. In fact, Bochner's formula is proved for Ricci 
fiat graphs by Chung- Yau [3], but it is hard to get the volume control or the 
Poincare inequality in that case. Lin- Yau |40] proved Bochner's formula for 
general graphs, but this does not lead to a version of Yau's gradient estimate 
that is nicely analogous to the case of Riemannian manifolds with nonnegative 
Ricci curvature. In the Cayley graph case, Kleiner [5H] obtained the Poincare 
inequality, but for non-abelian groups, Bochner's formula is unavailable (con- 
sider the free group case). For some special graphs which can be embedded 
into a surface with nonnegative sectional curvature in the sense of Alexandrov, 
Hua-Jost-Liu [26] proved the volume doubling property and the Poincare in- 
equality but Bochner's formula. It seems that Bochner's formula is sensitive to 
the local structure, but the volume growth property and the Poincare inequal- 
ity are not, c.f. [TS]. Hence the abelian groups are very suitable candidates 
for the application of Bochner's formula and Yau's gradient estimate. In ad- 
dition, Alexopoulos' theorem [3] is more general than ours, but it depends on 
the embedding of the nilpotent subgroup to simply connected Lie group. It 
seems hard to calculate the precise dimension of polynomial growth harmonic 
functions. We shall give the dimension estimate in an upcoming paper. 



dim HP k (R m ) 



2 Preliminaries and Notations 

Let G be an abelian group. It is called finitely generated if it has a fi- 
nite generating set. In this paper, we assume that any generating set S = 
{si , S2 , - ■ - ,S2i} of G is symmetric, i.e. S — —S, or more precisely s, = 
— Sj+;,1 < i < I, but we do allow that elements of S, are repeated, that is 
possibly Si — Sj for some i ^ j. We also allow G S. For any finitely gener- 
ated abelian group G with a generating set S, we have the associated Cayley 
graph (V, E) for which V = G, and xy £ E, (or x ~ y) if y — x 6 S, for 
x, y G V. The duplicity of elements in S produces multiedges between vertices, 
and G S makes self-loops. The vertices x and y are called neighbors if x ~ y. 
The degree of a vertex x is the number of its neighbors. Note that all vertices 
in (G, S) have the same degree (JS 1 . For the lattice Z n with the standard gener- 
ating set S a = {ei}f" l5 where = (0, • ■ • , 1, ■ ■ • , 0) is the i-th unit vector and 
ei = —e.i +n . 1 < i < n, we obtain the standard integer lattice in 1R™. This is the 
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object Heilbronn [22] and Nayar [JT] studied. In this paper, we consider the dis- 
crete harmonic functions on Cayley graphs of G with arbitrary finite generating 
set S. 

The Cayley graph of (G, S) is endowed with a natural metric, called the 
word metric (c.f. [6]). For any x, y E G, the distance between them is defined 
as the length of the shortest path connecting x and y, 

d (x, y) :— infjfc E IN | 3 x — xq ~ x\ ~ • • • ~ Xk = y}- 

Denote by Bjr(x) := {y E G d s (y,x) < r} the closed geodesic ball centered 
at x of radius r (r > 0). The volume of B^(x) is := )JG n B^(x), i.e. 

the number of vertices contained in B^'(x). For the subset SI C G, <i s (x, f2) := 
M{d s (x,y) y E £1} for any x E G, dfl := {z E G \ d s (z,fl) = 1}, and 
= QU dft. For any function / : 0, — > R, the discrete Laplacian operator is 
defined on ft as (x E fi) 

L s f(x) = ^2(f(y)-f(x)). 

The function / is called discrete harmonic (subharmonic) on f2 if L s f{x) = (> 
0), for any x E ft. In the lattice case (Z n , S°), the Laplacian operator is 

2n 

L s °f(x)=J^(f(x + e i )-f(x)). 
i=i 

The gradient of / at x E Q is defined as |V s /|(x) = y/J2y~x(f(y)-f( x )) 2 - We 
also need the partial difference operator (for s E S and x € f2) 

<5 s /(x) := /(x + a) - /(a). 

For Cayley graphs of (G, Si) and (G, S 2 ), it is known (c.f. [Ml US) that 

Gid Sl (^y) <d S *{x,y) <C 2 d s >(x,y), 

for any x, y E G, where Gi and G2 depend only on Si, S2. They are bi-Lipschitz 
equivalent in the metric point of view. Hence for any x E G, r > 0, 

l^(x)l<|£^)|;^(*)|. 

By the fundamental theorem of finitely generated abelian groups [13 H3]; an y 
finitely generated abelian group G is isomorphic to the direct sum Z' m (J) @\ = {L qi , 
where m, I € IN, = p" 1 for some prime number and some dj € IN. Hence 
there exists a generating set S* = {e±, • • • , e2 m , wi, ■ ■ • , W21} (e, = — ei+ m , ujj = 
-Wj+z, for 1 < i < m, 1 < j < I) such that G is identified with Z m ©j =1 Z g< , 
where {ei • ■ ■ e2 m } generates the torsion- free part Z m and {wi • • • W21} gener- 
ates the torsion part (B l i=1 1i qi . For (G, S°) it is easy to see that Gi(m, S°)r m < 
\Bf(x)\ < G 2 (m, S*°)r m , and |B| r °(x)| < G 3 (m, S°)\Bf (x)\, for any x E G and 
r > 1. Hence by the bi-Lipschitz equivalence, for the Cayley graph of (G, S) we 
have 

Ci(m,S°,S)r m < |fif(x)| < G 2 (m,S' ) S)r ro , 

\B§ r (x)\ <C 3 (m, S°, S)\B?(x)\. (2.1) 

The volume growth property (|2.1[) is called the volume doubling property. 

In the sequel, for simplicity we shall omit S from our notation, e.g., B r {x) :— 
B^ (x), if it does not cause any confusion. Also harmonic functions on G mean 
discrete harmonic functions. And constants C may change from line to line. 
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3 Yau's Gradient Estimate 

In this section, we prove Bochner's formula on finitely generated abelian 
groups and derive Yau's gradient estimate analogous to the case of Riemannian 
manifolds with nonnegative Ricci curvature. 

Kleiner [55] proved the Poincare inequality on Cayley graphs of finitely gen- 
erated (not necessarily abelian) groups. Let (G, S) be the Cayley graph of the 
finitely generated abelian group 

G-Z'»0©LA,- (3-1) 

By the volume doubling property (|2.1[) . we obtain the uniform Poincare inequal- 
ity. 

Lemma 3.1 (Poincare inequality, [IS]). Let (G, S) be the Cayley graph as &3.1\) . 
Then there exists a constant G(m, S) such that for any p G G, R > and any 
function f : B^^p) — > R we have 

E (f(x)-fB R ) 2 <CR 2 £ {f{x)-f{y))\ (3.2) 

xGB R (p) x,yeB 3R (p);x~y 

where f Ba = J2xeB R ( P ) f{x)- 

Note that by the independent works of Delmotte [17] and Holopainen-Soardi 
[21] the Moser iteration can be carried out for harmonic functions on graphs 
satisfying the volume doubling property and the Poincare inequality. The Cac- 
cioppoli inequality for harmonic functions was obtained for general graphs with 
bounded degree (c.f. [H M HQ)- 



Lemma 3.2 (Caccioppoli inequality). Let (G,S) be the Cayley graph as 13.1 
For any harmonic function f on Bqr(p) C G, R > 1, it holds that 



E iv/i 2 w<^ E A*). ( 3 - 3 ) 

xEBr(p) x£B 6R (p) 



where C = C(S). 



The Moser iteration implies the Harnack inequality for positive harmonic 
functions. 

Lemma 3.3 (Harnack inequality). Let (G,S) be the Cayley graph as 
Then there exist constants Ci(m, S) and C2{m,S) such that for any p £ G, 
R > 1 and any positive harmonic function f on Bc x r{p) we have 

max f < Ci min f. (3.4) 

B R (p) B R {p) 

The mean value inequality follows from one part of Moser iteration (c.f. 

[HOUHI])- 

Lemma 3.4 (Mean value inequality). Let (G,S) be the Cayley graph as 
Then there exists a constant C\(m, S) such that for any p G G, R > and any 
harmonic function f on Br{p) we have 

1 HU n x£B R (p) 



G 



The following Liouville theorem is a corollary of the Harnack inequality (|3.4p . 

Lemma 3.5 (Liouville theorem). Let (G, S) be the Cayley graph as i3.1\) . Then 
any nonnegative harmonic function f on G is constant. 

Bochner's formula has been obtained on Ricci flat graphs (c.f. [HI HO])- For 
the case of Cayley graphs of finitely generated abelian groups, we present the 
proof here for the convenience of readers. 

Lemma 3.6 (Bochner's formula). Let (G,S) be the Cayley graph as &3.1\) and 

f be a harmonic function defined on Bi(x), for i£G. Then 

L s \Vf\ 2 (x) > 0. (3.6) 

Proof. Let us denote S = {si, S2, ■ ■ ■ , S21}, where Sj = — Sj+j, 1 < i < I. Then 

2/ 

|v/| 2 (x) = J2(f(y) - /W) 2 = E l<WI 2 (*)- 

y~a; i—1 

Without loss of generality, it suffices to prove 

L s \S s J\ 2 (x)>0. 



l s \s s j\ 2 (x) = Ed^/I 2 ^)-!^/! 2 ^)) 

= ^2(f(y + s 1 )-f(y)f-2l\6 s jf(x) 

yr^x 

> ^[J2f(y + ^)-f(y)] 2 -mj\ 2 (^ (3.7) 

yr^x 

= ^E K x + s > + s i) - E /(y)] 2 - 2 ^i/l 2 0*0 

i—1 y^x 

= ~[ E /W-E^f)] 3 -^/!^) 

= ^[2;/(x + Sl )-2//(x)] 2 -2^ Sl /| 2 (x) (3.8) 
= 0, 

where we use the Holder inequality in (|3.7[) and the harmonicity of / in (3.8). □ 

Combining Bochner's formula with previous results, we obtain Yau's gradient 
estimate for positive harmonic functions. 

Proof of Theorem ll.il We choose C\ — 7C[ , where C[ > 1 is the constant 
Ci(m,S) in Lemma [3.31 Bochner's formula (|3.6[) implies that |V/| 2 is a sub- 
harmonic function on Bc 1 r{x). Then the theorem follows from the mean value 
inequality (|3.5|) (since 5 Si f is harmonic for any £ S), the Caccioppoli in- 
equality (|3.3D . the volume doubling property (|2.1j) and the Harnack inequality 
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(EH , 



- R 2 \B R {x)\ 5- /2(2/) 



□ 



Corollary 3.7. Le£ (G, S 1 ) &e i/ie Cayley graph as VS. 1}) . There exist constants 
Ci(m, S) and Ci{m, S) such that for any x £ G, R > 1 any harmonic function 
f on Bc-ir(x) we have 

|V/|(aO<^ osc : /, (3.9) 

ti B Ci r{x) 

where osc Bc , ir(x) f = max BciJi(l) / - mm BaiB ( x ) /. 

Proof. It suffices to choose /— ming c fl ( x ) /+£ (small e) as the positive harmonic 
function in Theorem ll.il □ 



4 Polynomial Growth Harmonic Functions are 
Polynomials 

The fundamental theorem of finitely generated abelian groups implies that 
for any finitely generated abelian group G there exists a generating set 5° = 
{ei,-'- ,e 2 m,wi,--' ,w 2 i} (e; = -e i+m , wj = -w J+l , for 1 < % < m, 1 < j < I) 
such that G = Gi © G 2 is identified with Z m 0®- =1 Z gi , where {ei ••• e 2m } 
generates the torsion- free part G\ = Z m and {wi ■ ■ ■ W21 } generates the torsion 
part G 2 = ©i =1 Zg 4 . For some fixed p G G, we denote by # d (G, 5) = {it : G ->• 
R I L s it = 0, 3G s.i. |tt(jc)| < C(d s (p,x) + l) d } the space of polynomial growth 
harmonic functions on G with growth rate less than or equal to d. For fixed S°, 
Gi is identified with Z m , then we denote by P d {Z m ) := P d {R m ) the space of 
polynomials in R m restricted to the lattice Z7 n with degree less than or equal 
to d. 

First, we consider the easy case that G is torsion-free, i.e. G = Z m . The 
following theorem generalizes the classical theorem of Heilbronn [22] . 

Theorem 4.1. Let (G, 5) be the Cayley graph of a torsion-free finitely generated 
abelian group, G = Z m . TTiera polynomial growth harmonic functions on G are 
polynomials, i.e. 

H d (G, S) C P d (Z m ). 

Proof. Let S 10 = {ei, • • • , e 2m \ be the standard basis for Z m , S — {si, • • • , S2/} 
be the generating set for the Cayley graph (G, S). Then 

21 

ei = ^a^Sk, 

k=l 
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where a$ € Z for 1 < i < 2m, 1 < k < 21. Let 

6if(x) :=5 e J{x) = f{x + e i )-f{x), 
for 1 < i < 2?n. Since G is abelian, it is easy to show that 

L S SJ = 5 t L s f = 0, 

for the harmonic function /. 

We claim that SJ G H d ~ 1 (G, S) if / £ ff d (G, 5). Although for any a; G G, a; 
and x + may not be neighbors in the Cayley graph (G, S), there exists a path 
from x to x + ei = x + Y^k=i a i s k, he. x = xq ~ ~ ■ ■ • ~ x t = a; + e*, whose 
length is t = J^lli \a k i I < C(S°, 5). Note that / G H d (G, S) <^=> osc Br(p) / < 
Gi? d for some fixed p G G and sufficiently large R > Rq. For any x G G, 
1 < i < 2m, 

|<fc/IO»0 = |/(x + e0-/(x)| 

< l/(x + e,) - /(x t _i)| + |/(x t _r) - /(x t _ 2 )| + ■ ■ ■ + |/(xr) - /(x)| 

t-i 

< — osc J, 

R B 2 r(x) 

for R > Ri(S°,S), since Xj G Bc(x), for < J < i — 1. The last inequality 
follows from Corollarv l3.7l For any x G Br(p), we have B 2 r(x) C B 3 r(p). Then 

G 

osc < 2 max l^/lfx) < — osc / < Gi? d_1 , 

for R> Ri. This proves the claim. 

Hence by taking finitely many times partial differences, we obtain 

for any k\ + k 2 H 1- fem > M + 1 an d / G H d (G, S), where [d] is the maximal 

integer not exceeding d. By the basic difference equation theory or Lemma 2.13 
in Nayar [IT], we conclude that / is a polynomial. □ 

Then we can prove the generalized Heilbronn's theorem. 

Proof of Theorem 1.2. It suffices to show that / is constant on G 2 , i.e. f(x + 
w) = /(x), for any x G G, w G G 2 . 

Let 5° = {ei,--- ,e 2m ,wi,-- - ,w 2 i} (ej = -e i+m , Wj = -Wj+i, for 1 < 
i < m, 1 < j < I) such that G = Gi © G 2 is identified with Z m 0e^ =1 Z (?i , 
where {ei---e2 m } generates the torsion-free part G\ = Z m and {w,, Wj+i} 
generates Z g4 . Let S w .f(x) := f{x + — f{x). The same argument as in 
the proof Theorem oWlies that 8 W] f G H d - 1 {G,S) if / G H d (G,S). Then 
(5^ + V = 0, for / G H d {G, S). For fixed x G G, lifting Z 9 . to Z, we obtain that 
/ is a polynomial on Z. Since /(x + (r + kqj)wj) = f(x + rwj) for any fc, r G Z, 
the periodic polynomial / must be constant, i.e. /(x + ruty) = /(x), for any 
r G Z. Since this is true for any 1 < j < i, we obtain /(x + w) = /(x), for any 
x G G, u; G G 2 . Then it is easy to see that L" G i s /(xi) = for any xi G Gi if 
L s f{x) = for any x G G. Hence H d (G,S) = H d (G 1 ,Tt Gl S) C pW(Z m ). □ 
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By the Harnack inequality, we reprove Nayar's theorem. 

Proof of Theorem 1.3. It suffices to show HM d (G, S) C H d (G, S). Without loss 
of generality, for / G HM d (G, S) satisfying f(x) > —C(d(p, x) + l) d , we need 
to prove that f(x) < C(d(p,x) + l) d , for some C. For simplicity, we assume 
f(p) — 0. Let G\ be the constant for the Harnack inequality in the Lemma 13731 
Then for any x 6 Br(p), R > 0, it is easy to see that Bc 1 r(x) C B^c 1 +i)b^{p)- 
Moreover 

f(y) > -C(d( P , y) + l) d > + 1)R + l) d > -CR d , 

for y e B (Cl+l)R {p), R>1. That is f(y) + CR d > on B Ci r(x). The Harnack 
inequality (|3.4|l implies that 

f(x) + CR d < C(f( P ) + CR d ) = C 2 R d - 

Then we have 

f(x) < CR d , 

for x £ Br(p), R > 1. Hence there exists a constant C such that /(x) < 
C(d(p, x) + l) d . □ 



5 Calculating the Dimension 

For calculating the dimension, by Theorem 11.21 it suffices to consider har- 
monic polynomials on a torsion-free finitely generated abelian group. Let (G, S) 
be the Cayley graph of G = Z m . There exists a generating set S° — {e\, • • • , e2«} 
such that G is identified with Z m in R m . For k e IN U {0}, we denote by 
P fe (R™ 1 ) the space of polynomials on W 71 of degree less than or equal to k, by 
Pi := P k (Z m ) := {u : Z m ^ R | = G P fc (R m )} the space of the 

restriction of polynomials on R m to Z m of degree less than or equal to k. We 
denote by := HP k (R m ) the space of harmonic polynomials on R™ 1 (Au = 0) 
of degree less than or equal to k. For the Cayley graph (G, S) which is identified 
with (Z m ,S), we set D k Sm := HP k {Z m ,S) := {u E P k {Z m ) \ L s u = 0}. In 
order to calculate the dimension of discrete harmonic polynomials, we make the 
dimension comparison between D k s and R k n . It is well known for harmonic 
polynomials on R m (see [33]) that 

(k > 1, dimi?° n = 1) 

(*>1) 

(k > 0) 

(k > 2) 
Lemma 5.1. 

dim£»| 1 = dim^ = 2, (5.1) 

for any S, k > 1. 



dimR k n 
dimi?^ 
dimP^ 
dim P k 



m + k — 1 

k 



m + k — 2 
jfe-1 



dimR k n _ 1 +dimP^ 1 , 



= dimP^R" 1 ) =^ 



m + i — 1 



dim P!; + dim P. 



fc-2 
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Proof. Let S — {o»}fLn en G Z, and aj = — for 1 < i < L Since S is a 
generating set, at least one at ^ 0. For any / g Z)| \, f = bkX k + bk-\x k ~ 1 + 
• • • + bo, where bo, b\, ■ ■ ■ ,bk £E. By Taylor expansion, the difference equation 
L s f = Q reads 

n=l i=l v ' 

Comparing the degree of polynomials, we have 

{b k x k y = o, 

for x £ 1i. Hence, k < 1, i.e. / is linear. □ 
Lemma 5.2. 

dim^ m <dimi4, (5.2) 

for any S, k > 0, m > 1. 

Proof. It is easy to see that -Djg m = = const. We apply an induction 
argument on fc. It suffices to prove (|5.2p for /c = Z if (|5.2D is true for all fc < I — 1. 
We pick ei g S* and define Sif(x) = f(x + ei) — f(x), for any function / on G. 
Since G is abelian, 8\L S = L Si, then 6i is a well defined linear operator, 

5i:Dl m ^D k s ;i. 

By linear algebra, 

dimDg m = dimker<5i + dimim Si, (5-3) 

where keriJi and im <5i are Kernel and Image of Si. There is a natural projec- 
tion P : Z m -> Z" 1 - 1 , for any a; = (xi,x 2 ,--- ,z m ) <E Z m , Px = x', where 
x' = (x2, • ■ ■ ,x m ). For any generating set S for Z7 n , S" = {s' \ s E S} is the 
generating set of Z m_1 . Moreover, S^' :— P l S is the generating set for Z m_t , 
1 < t < m- 1. 

For any / € ker<5i, i.e. <5i/ = 0, then f(xi, x 2 , ■ • • ,x m ) = g(x 2 ,--- ,x m ). 
Hence 

= L s /(x) = ]T(/(x + ,s)-/(x)) 

= ^( 9 (,' + s ')- s ( I '))^ s s(4 

That is ker^i = Z5|, m _!- Hence (|5-3[) implies that 

dimP| m < dimi^l,^.! + dimD*^, 



11 



for any S, k > 1, m > 1. Then it follows that 

dim^ m < dimZ^m-i+dimDijl 



< dim + dim D^7 1 ro _ 1 + dim 

< 

m 

< dim D^(m-i) i + dim D s7m- 4 ) i 



m 

< 2 + ^dimi?^ 1 

m 

= dim + ^ dim flj 

i=2 

= dim R.L, , 



where we use Lemma 15.11 the inductive assumption (15.21) for k < I — 1 and some 
facts in ft™. 

□ 



Now we can prove the main Theorem ll.4l 

Proof of Theorem \l-4\ It suffices to show that 

dimL>|, m = dimi?* , 

for any S, k > and m > 1. We may assume k > 2, otherwise it is trivial. By 
S = —S, the Laplacian operator L s is a linear operator, 



tS . pk , pk — 2 

J • ill ' J> 1 • 



Since kerL s ~ Z?| m , we have 

dim P r k n = dim kcr L s + dim im L s 

< dimZ?| m + dimP„ fe - 2 

< dim R k n + dim P^~ 2 
= dim P k 

which follows from (15.21) . Hence dim_D| m = dimi?^. □ 

Remark 5.3. In the above theorems, all the inequalities for dimension com- 
parison are equalities. Hence, we obtain that 8\ and L s are surjective linear 
operators. In fact, Heilbronn |22] used the technical lemma in difference equa- 
tion theory (see [MIDI]) that L s : P^ -> P^ 2 is surjective to calculate the 
dimension. Conversely, by dimension comparison, we obtain a more general 
difference equation lemma. 



Corollary 5.4. Let (Z m ,S) be the Cayley graph ofZ m . Then 

Si ■ D k S m D k s -^, 



and 

r S . pk . pk—2 

• m m 

are surjective linear operators. 
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